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NOMENCLATURE 



English 


Letter Symbols 


A 


Area (ft^) 


A1 


- Estimated Area for Boiling in First Superheater 
Pass (ft^) 




2 

Frontal Area Blocked by Tubes and Fins (ft ) 




2 

Bare Tube Area (ft ) 




2 

“ Heat Exchanger Frontal Area (ft ) 


^fin 


2 

- Fin Area (ft ) 


A^i^ 

1 


2 

Cross-Sectional Area for Fluid Flow (ft ) 

Inside Area Required for Boiling in First 

2 

Superheater Pass (ft ) 


A^^ 

1 


Inside Area Required for Superheating in First 

2 

Superheater Pass (ft ) 


A. 

IP 


2 

Inside Heat Transfer Area Per Pass (ft ) 


A . 
min 


2 

- Minimum Cross-Sectional Area for Gas Flow (ft ) 


^op 


2 

- Outside Heat Transfer Area Per Pass (ft ) 




2 

Total Heat Exchanger Inside Area (ft ) 


to 


2 

- Total Heat Exchanger Outside Area (ft ) 


C 


Constant 


*“bhp 


Correction Factor to Gas Turbine BHP for Duct Loss 


C 

max 


- Maximum Heat Capacity (Btu/hr-F) 


C . 
min 


- Minimum Heat Capacity (Btu/hr-F) 


C 

Pf 


- Specific Heat of Water/Steam (Btu/Lfcin-F) 


C 

pg 


- Specific Heat of Gas (Btu/lbm-F) 


*“sfc 


Correction Factor to Gas Turbine BHP for Duct Loss 



11 



I 



- Fin Outside Diameter (ft) 

d^j^ - Diameter of Fin Base (ft) 

dj^ - Inside Tube Diameter (ft) 

d^ - Outside Tube Diameter (ft) 

d^ - Fin Root Diameter (ft) 

f - Friction Factor 

2 

G - Gas Flow Rate Per Square Foot (Ibm/hr-ft ) 

2 

” Maximum Gas Flow Rate Per Square Foot (Ibm/hr-ft ) 

- Enthalpy of Steam at Turbine Inlet (Btu/lbm) 

- Turbine Exhaust Steam Enthalpy (Btu/lbm) 

h2g “ Turbine Exhaust Steam Enthalpy Assuming Isen- 

tropic Expansion (Btu/lbm) 

h^ - Enthalpy of Saturated Water (Btu/lbm) 

hgj^ - Enthalpy of Water at Heat Exchanger Inlet (Btu/lbm) 

h-„ - Enthalpy of Water at Heating Section Outlet 

^ (Btu/lbm) 

h-_ - Enthalpy of Steam at Boiling Section Outlet 

(Btu/lbm) 

h^^ - Enthalpy of Steam at Superheater Outlet (Btu/lbm) 

h^g - Enthalpy Increment for Boiling (Btu/lbm) 

h^^ - Enthalpy of Heat Exchanger Feedwater (Btu/lbm) 

hp - Heat Exchanger Inside Heat Transfer Coefficient 

^ 2 

in Heating or Superheating Section (Btu/hr-ft -F) 

h^p^ - Heat Exchanger Inside Heat Transer Coefficient 

in the Two-Phase Region (Btu/hr-f t^-F) 

j - Heat Transfer Colburn j -factor 

k - Thermal Conductivity of Gas (Btu-hr-f t-F) 

g 

k^ - Thermal Conductivity of Steam/Water (Btu/hr-f t-F) 
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N,NTU 

N 

s 

N,. / 
t/r 

P 

^gt 

"st 

^tot 

Q 

Qb 

Qh 

S 

"^sh 



- Thermal Conductivity of Heat Exchanger Tube 
Wall (Btu/hr-ft-F) 

Fin Height (ft) 

Length of Cut from Fin Tip (ft) 

Tube Length (ft) 

- Steam/Water Flow Rate (Ibm/hr) 

- Gas Flow Rate (Ibm/hr) 

Gas Flow Rate in Boiling Section of First 
Superheater Pass (Ibm/hr) 

- Gas Flow Rate in Superheating Section of First 
Superheater Pass (Ibm/hr) 

Number of Passes 

Number of Transfer Units 

Number of Fins Per Inch 

Number of Segments in One Fin 

Number of Tubes Per Row 

Pressure (psia) 

Steam Pressure at Steam Turbine Inlet (psia) 
Steam/Water Pressure in Heat Exchanger (psia) 

- Gas Turbine Horsepower 

- Steam Turbine Horsepower 
Total System Horsepower 

- Heat Transfer Rate (Btu/hr) 

- Heat Transfer Rate in Boiling Section (Btu/hr) 

- Heat Transfer Rate in Heating Section (Btu/hr) 

- Heat Transfer in a Heat Exchanger Pass (Btu/hr) 

- Heat Transfer Rate in Superheating Section (Btu/hr) 
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Qrb 

q" 

R 

o 

s 

"fg 



^fb 

Tfi 

^f2 

^f3 

Tf4 

4 

T n 

gi 

T o 
g2 

^g3 

^gb 



Required Heat Transfer Rate in Boiling Section 
of First Superheater Pass (Btu/hr) 

Heat Flux (Btu/hr-ft^) 

2 

Thermal Resistance (hr-ft -F/Btu) 

2 

- Heat Exchanger Outside Resistance (hr-ft -F/Btu) 
Entropy (Btu/lbm-F) 

Entropy of Steam at Turbine Inlet (Btu/lbm-F) 
Entropy of Saturated Water (Btu/lbm-F) 

Entropy Increment for Evaporation (Btu/lbm-F) 

Tube Spacing Normal to Gas Flow (ft) 

- Tube Spacing Parallel to Gas Flow (ft) 

Steam Temperature at Turbine Inlet (F) 

Steam/Water Bulk Temperature (F) 

Water Temperature at Heat Exchanger Inlet (F) 

Water Temperature at Heating Section Outlet (F) 
Steam Temperature at Boiling Section Outlet (F) 

- Steam Temperature at Superheating Section 
Outlet (F) 

Fin Thickness (ft) 

Average Gas Temperature at Heat Exchanger Inlet (F) 

- Average Gas Temperature at Boiling Section 
Inlet (F) 

- Average Gas Temperature at Heating Section 
Inlet (F) 

- Average Gas Temperature at Heat Exchanger 
Outlet (F) 

Gas Bulk Temperature (F) 
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T ^ 

gf 

■^sat f 
T 

wo 

U . 
oi 

w 

s 

Wt 



X 

a 

^2 

^2s 



Gas Side Film Temperature (F) 

Temperature of Saturated Water (F) 

Average Outside Tube Wall Temperature (F) 

2 

Overall Heat Transfer Coefficient (Btu/hr-ft -F) 
Fin Segment Width (ft) 

Steam Turbine Work (Btu/lbm) 

Isentropic Steam Turbine Work (Btu/lbm) 

Steam Quality 
Average Steam Quality 
Turbine Exhaust Steam Quality 

Turbine Exhaust Steam Quality Assuming Isentropic 
Expansion 



Dimensionless Groups 

N - Nusselt Number 

u 

P - Prandtl Number 



r 

- Reynolds Number 
S^ - Stanton Number 



Greek Letter Symbols 



AP 

AT 

£ 

e. 



oa 



Pressure Change (psia) 

Temperature Change (F) 

Effectiveness 

Effectiveness of Boiling Section of First 
Superheater Pass 

Effectiveness Overall for a Heat Exchanger Section 
Effectiveness for a Heat Exchanger Pass 
Viscosity at Bulk Temperature (Ibm/ft-hr) 
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’'st 

p 

Pv 



- Viscosity at Tube Wall Temperature (Ibm/ft-hr) 

- Fin Efficiency 

System Thermal Efficiency 

- Steam Turbine Efficiency 

3 

Density (Ibm/ft ) 

Density of Saturated Water (Ibm/ft^) 

3 

Density of Saturated Vapor (Ibm/ft ) 
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INTRODUCTION 



I . 



A . BACKGROUND 

The U. S. Navy has made a commitment to gas turbine 
propulsion systems, with two classes of gas turbine ships 
entering the fleet (DD-963 and FFG-7) and one class pro- 
posed (DDG-47) at the time of this writing. The DD-963 class 
has two shafts with two General Electric LM 2500 gas tur- 
bines powering each shaft. The FFG-7 class is a single- 
shaft ship with two LM 2500 gas turbines. The DDG-47 class 
propulsion plant, as proposed, will be similar to that of 
the DD-963 class. 

Currently , means for conserving fuel are being sought 
in all sectors of the economy and government. In most of 
our steam turbine-powered ships fuel conservation must be 
pursued largely through operating practices such as reduced 
steaming time or reduced speed when steaming. Several 
alternatives exist, however, for fuel consiimption reduction 
for the gas turbine propulsion system. These alternatives 
involve combining another type of propulsion system with 
the gas turbine in order to reduce the overall fuel consximp- 
tion by either operating the combined systems in parallel 
or by operating the less expensive system in the cruise 
mode. Some examples of these combinations are: 

Gas Turbine and Diesel (CODOG, GODAG) 

Gas Turbine and Steam (COGAS) (steam system not 
separately fired) 
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Gas Turbine and Steam (COSAG) (steam system 
separately fired) 

The combined Gas turbine And Steam turbine (COGAS) system 
is probably the best of these alternatives with respect to 
fuel consumption, initial cost, maintenance, and manning 
requirements . 

The COGAS ship propulsion plant uses the exhaust heat 
from a gas turbine to generate steam in a waste heat re- 
covery unit (WHRU) . This steam is used to drive a steam 
turbine which is operated in parallel with the gas turbine, 
for ship propulsion. A common reduction gear connects both 
turbines to the propeller shaft. Figure 1 illustrates the 
COGAS system. 

Several studies have examined the potential of the COGAS 
system for improving propulsion plant fuel efficiency. 

Giblon and Rolih [Ref. 1] studied the commercial ship 
application of COGAS. References 2 and 3 investigated the 
feasibility of the COGAS system for military ship application. 
In all three studies the COGAS system was designed around 
a gas turbine brake horsepower of about 23,000. Additionally, 
all three employed a drum-type boiler as the waste heat 
recovery unit. 

In this study, COGAS system designs were considered for 
a destroyer-sized ship with two shafts and two LM-2500 gas 
turbines for each shaft. One gas turbine was equipped with 
a WHRU. In the COGAS mode of operation the WHRU-equipped 
gas turbine was operated in parallel with the steam turbine. 
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CONDENSATE PUMF^ 



powering one shaft. The other shaft was allowed to drag. 

It was assumed that the principal employment of the COGAS 
system would be for "cruise" conditions. That is, the COGAS 
system would be operated primarily for steady steaming. The 
horsepower range over which the COGAS system could be 
operated is, theoretically, limited only to the range of 
powers over which the gas turbine can be operated. Thus, 
the power range would vary from the COGAS system power 
output with the gas turbine at idle to the system power 
output attained with the gas turbine operating at its maxi- 
mum continuous power rating. 

When the COGAS system operating range is selected, the 
desired goal of fuel conservation must be balanced against 
system reliability and routine maintenance requirements . 

If only one gas turbine is equipped for COGAS operations, 
the system downtime will increase as the turbine is operated 
at higher power settings. In order to increase the time 
between major gas turbine maintenance requirements, two 
COGAS systems could be installed in a twin-shaft ship. The 
increased overall reliability which this redundancy would 
afford, however, would be at least partially offset by the 
higher initial cost, the increased manning requirements, 
and the added weight of two systems. 

A once-through counter-cross-flow heat exchanger was 
selected as the WHRU. Reference 6 states some of the 
principal advantages of the once-through steam generator 
over the drum-type boiler. These advantages are as follows: 
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lower initial cost, faster response, more compactness, 
and lower operating cost. Additioitally , the once-through 
unit provides a somewhat simpler model than the drum- type 
heat exchanger. A schematic representation of a once- 
through WHRU installed in the exhaust ducting of a gas 
turbine is shown in figure 2. 

B . OBJECTIVES 

There were five main objectives for this thesis. 

1. Develop a design model for the once-through WHRU. This 
model was to size a WHRU and predict is performance given 
initial conditions such as feedwater inlet temperature, 
steam outlet temperature, WHRU operating pressure, inlet 
and outlet gas conditions, and gas flow rate. 

2. Consider several sets of initial conditions and design 
constraints (e.g., maximxam frontal dimensions) and produce 
WHRU designs for each set. 

3. Evaluate these designs and integrate them with a simple 
COGAS system model in order to test the design model and 

to develop a clearer understanding of the WHRU performance 
in the COGAS system. 

4. From the evaluated designs, predict the performance 
of the COGAS system. 

5. Provide a framework for future studies involving a 
detailed, complete COGAS system cycle analysis and 
optimization . 
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OMCE-THROUGIl V/ASTE MEAT RECOVERY UNIT 



II. MODEL DESCRIPTION 



A. OVERVIEW AND INITIAL CALCULATIONS 

The waste heat recovery unit (WHRU) is modeled as a 
once-through counter-cross-flow heat exchanger. The model 
is applied in a computer program for use in heat exchanger 
design for a combined gas and steam (COGAS) turbine pro- 
pulsion plant system. The same model could be adapted 
for use in a complete cycle analysis of the COGAS system. 

The WHRU model is divided into three principal sections: 
heating, boiling and superheating. The specified initial 
conditions for the model are: (1) gas temperature into 

the heat exchanger , (2) minimum average gas exit 

temperature from the heat exchanger (Tg^) , (3) gas flow 

rate (*^1 , (4) water inlet temperature (T^^) , (5) water/ 

Steam pressure (P^) • Useful background material for the 
model formulation was obtained from references 1-5. 

Two simplifications should be noted at this point. First, 
the water-side pressure drop in the WHRU is neglected. 

Second, gas temperatures, evaluated at several points in 
the model, are average temperatures. Specifically, the 
average inlet and outlet gas temperatures for the heating 
and superheating sections and the average inlet and outlet 
gas temperatures for each pass in the boiling section are 
evaluated. In each of these calculations the gas tempera- 
ture distribution is neglected. Clearly, even if the gas 
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turbine exhaust flow is uniform in temperature as it 
approaches the WHRU, the gas temperature will not be uni- 
form after the first WHRU pass is encountered. Due to the 
varying gas-to-fluid temperature difference across each 
pass the gas temperature will be distributed non-uniformly 
after the first WHRU pass. One can, however, calculate the 
average gas temperature approaching a particular pass 
based on the heat release from the gas in the previous pass 
and the average gas temperature entering the previous pass . 
Averaged gas temperatures calculated in this manner are 
used throughout the model. 

Prior to the heating section, some initial calculations 
are performed to determine the water/steam mass flow rate 
and to set the pinch point. In order to calculate the 
water/steam mass flow rate an energy balance is performed, 
using the starting conditions, as follows; 



Q = A C (T 1 -T 4 ) 
g pg gl g4 



where C is evaluated at the average gas temperature 
PU 

T = (T , +T .)/2. The fluid mass flow rate is 

gavg gl g4 

calculated from 



m 



f 



Q 

^f4 "^fl 



where h^j^ and h^^ are inlet and outlet fluid enthalpy 



24 



respectively. Intermediate fluid and gas temperatures are 
calculated as follows: 

1. The heating section is assumed to heat the water from 
inlet conditions to the saturation temperature, ^ 

heat balance is performed on the water side yielding a 
heating section heat transfer of 



Qh = 



The required average gas temperature entering the heating 
section is then calculated from 



T , = T . + 

g3 g4 



m C 

g pg 



2. In the boiling section it is assumed that sufficient 
heat is added to produce saturated vapor at the outlet. 

An energy balance on the steam- side provides the total heat 
transfer in the boiling section. 



Qb m^(h^2"hf2^ 



where h^^ the enthalpy of saturated vapor at the 
pressure specified. The required average gas temperature 
entering the boiling section is 



T - = T , + 

g2 g3 



ih C 

g pg 
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3. In the superheating section, the steam temperature is 
raised from the saturation temperature to the outlet 
temperature specified in the initial conditions. The heat 
transfer in the superheating sections is 



°sh ' 



where h^^ is the enthalpy of steam corresponding to the 
steam outlet temperature The gas inlet temperature 

T , is known from the initial conditions. 

gi 

The "pinch point", for the purposes of this model, is 
defined as the difference between the gas temperature enter- 
ing the heating section and the water temperature leaving the 
heating section (saturation temperature) . This temperature 
difference will normally correspond to the smallest tempera- 
ture difference between gas and fluid in the heat exchanger. 
The only other likely location for the smallest temperature 
difference between gas and fluid is at the gas inlet, or 
superheater outlet. Figure 3 depicts the typical gas and 
fluid temperature distribution in the WHRU. 

The model must provide some internal control over the 
pinch point temperature difference as defined above. This 
is desirable because the gas-fluid temperature difference 
at the fluid saturation point is indicative of the amount of 
heat transfer area which will be necessary to heat the in- 
coming water to the saturation point. As the temperature 
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TYPICAL V/HRU GAS-FLUID TEMPERATURE DISTRIBUTION 




PERCENT OF TOTAL \mR\J LENGTH 



FIGURE 3 
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difference at the pinch point increases the heating portion 
of the heat exchanger becomes more efficient from the point 
of view of heat transfer area required. The minimum pinch 
point temperature difference is set at 25 F. After the 
interim temperatures are established the pinch point temp- 
erature difference is checked. If the difference is less 
than 25 F, the gas temperature at the heating section inlet 
is reset to 



To = T^_ + 25. 

g3 f 2 



From the overall energy balance 



C m (T , -T •) = m^(h_. - h_, ) 

pg g gl g4 f f4 fl 



and the energy balance across the heating section 



C m(To“T.) = m^(h_o“hp,), 

pg g g3 g4 f f2 fl 



a revised gas outlet temperature may now be calculated from 



T o ” ctT , 

T = -2l 

^g4 1-a 



where 



*'f2 - •'fl 
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After this new gas inlet temperature is established, all 
initial calculations are performed again, to establish a 
new water/steam mass flow rate and to recheck the pinch 
point temperature difference. Once these calculations are 
complete, the geometric parameters for the heat exchanger 
must be specified. 

B. GEOMETRY 

A fin-tube with helically curved extended surfaces on 
circular tubes is assumed for this model. The fins are 
segmented. The tubes are configured in banks of one row 
each and the rows are staggered. The base tube surface is 
taken from Ref. 7. The description of the finned tubes 
is as follows. 

d^ = tube inside diameter = 1.86 in. 

d = tube outside diameter = 2.00 in. 

o 

d = fin root diameter = 2.00 in. 

r 

Ng = fins per inch = 5.94 

a = fin height = 1.015 in. 

= length of cut from fin tip = 0.82 in. 

d^ = fin outside diameter = 4.03 in. 

t^ = fin thickness = 0.048 in. 

w^ = fin segment width = 0.17 in. 

Ng = number of segments in 360 degrees = 38 

The finned tube configuration is shown in figure 4. The 
tube length, L, and number of tubes per row are 

chosen by the designer according to the frontal area, A^, 
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Figure 4; Segmented Fin-Tube Configuration 

desired. The tube layout is shown below in figure 5. 

The center-to-center tube spacing in the transverse direc- 
tion is 4.50 inches. The equilateral, staggered tube 
arrangement used in this model leads to a spacing normal 
to the gas flow, S^, of 4.50 inches and a spacing parallel 

to the gas flow, S , of 3.90 inches as shown in figure 5. 

P 

The heat exchanger height is established by the number 
of WHRU passes (rows) and the tube spacing parallel to the 
gas flow (Sp) • 
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4.50 in. 0.50 in. 



Figure 5 : Model Tube Layout 

In order to establish the minimum gas flow cross- 
sectional area the total "blocked" frontal area, A^, of 
the heat exchanger must be calculated from 

*b = '"t/r L L Nj L 1 tj , 

and the minimum gas flow area is 

A . = A. - A, 

min f D 

The total inside area available for heat transfer per pass 
is 




7T d . 
1 



L N 



t/r • 
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The gas side surface area available per pass for heat 
transfer is the sum of the fin surface area and the bare 
tube surface area per pass. The fin surface area per tube 
is calculated from 

^fin 

where d^j^ = d^ - The bare tube area per tube is 

A, . = TTd L-ird t^N^L. 

Dt o off 

Therefore, the total outside area per pass available for 
heat transfer is 



= N. L[N^(2 I + 2 t . + w t^) 

f s cs fc sf 



TT , ,2 j2. , 



A = , (A^. + A, 

op t/p fin bt 



The cross-sectional area for fluid flow is calculated from 



A 



ff 




N 



t/p • 



With the mass flow rates, terminal temperatures and heat 
exchanger geometry established, the remainder of the model 
may be solved for number of passes, actual interim tempera- 
tures, location of the fluid phase changes, and the gas 
side pressure drop. 
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C. HEATING SECTION 



The gas-side Reynolds number for the heating section is 
calculated initially using the gas bulk temperature to find 
the gas properties. With this Reynolds number, a j -factor 
and a friction factor, f, are obtained from polynomials fit 
to the data for tube layout number 5 in Ref. 7 (fig. 5) . 

The friction factor will be used later to calculate the 
gas side pressure drop in the heating section. The j -factor 
is related to the heat transfer coefficient h by the 

g 

following relationship. 




By introducing 



the previous expression can be written as 




N 

u 



J 




and 



N 



u 




where 
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N 



u 



h d 

_2 

K 



o 



Therefore, a relationship may be written for the heat 
transfer coefficient as follows: 



= 3 



K 

d 

o 



R 

eg rg 



The water-side heat transfer coefficient is calculated by 
using the Dittus-Boelter correlation 



h 



f 



0.023 





P 



0.4 

rf 



where all properties are obtained at the bulk temperature 
of the water in the heating section. Using the tube wall 
resistance together with h and h^, the overall heat transfer 
coefficient for the heating section may be written in terms 
of the inside area as 



U . 
oi 



In (d /d . ) 
1 + 2_i_ 

5T 2 7T K , L 

' f w t/p 



A, 

+ 



A n . h 
op 't g 



where 



n. = 1 - (1 - n^) ^fin . 

A 

op 
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The fin efficiency, n£, is calculated from the expression 

_ tanh ML 
'^f ML ' 



where 



M = /h P/KA 

g 

and, if the fin is approximated by a set of rectangular 

df - d 

strips extending from the tube wall, L = ^ . Now, 

the cross-sectional area. A, of a rectangular strip may 
be written as 



A = 




/ 



and the perimeter is 



P = 2 (w + t^) . 

s f 

So, using known geometric parameters and the thermal conduc- 
tivity of the fin metal, the parameter ML may be expressed 
as 



ML = C /h 

g 



where 
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c 



2 





/2(w^ + tj) 

V ”s 



NOW, 



= [1 - (1 - 



tanh C /h~ . 

2 _) -iin.] 

C /h~ ^op 

g ^ 



The number of passes required in the heating section is 
calculated, using the ef fectiveness-NTU method, in the 
following way. 

1. An average pass effectiveness for the heating section is 
calculated from the following expression for cross-flow 
effectiveness with both fluids unmixed obtained from 
Ref. 8, 



where 



, rSxp (-NCn) - 1, 

= 1 - exp[ - £ - ^ 3 



c = 



c . 

min 

C 

max 



C _ m^ 

p£ f 
C m ' 

pg g 



N = NTU = 



U . A. 
oi ip 

c . 

min 



n = N 



- 0.22 



Water and gas properties are taken at bulk temperatures. 
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2. An overall heating section effectiveness is calculated 
from 



£ 



oa 



1 - £ C . /C 
, p min max^n 

^ 1 - £ ^ 

E 

1 - £ C . /C 
. p min max.n 

^ 1 - £ ’ 



C . 
min 

C 



max 



where n = number of passes. This expression was obtained 
from Ref. 9. The overall effectiveness formulation was 
actually derived with the condition that fluids are mixed 
between passes. It can be shown [Ref. 9] , however, that 
the error is not large when the expression is used for the 
case where fluids are unmixed between passes. 

3. The expression for £^^ is now used in an iterative 
fashion starting with n = number of passes = 1. Each time 
£^^ is solved, the gas temperature into the heating section, 
T T , is found from 



T 



g3 



£ C T — C T 
oa min gl max g4 

£ C - C 
oa min max 



which is derived from 



£ 



oa 



C 

max 



C . 
min 






4. Once the inlet gas temperature is established for a 
particular number of passes, an energy balance may be 
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performed to solve for the total heat transfer which would 
take place in the heating section for that number of passes. 



= m C (T - - T J 
g pg g3 g4' 



5. The enthalpy of the water at the outlet of the heating 
section for a particular iteration may be found from 



h 



fL 




/ 



and this yields the temperature of the water at the outlet, 
Tf 2 » for the number of passes under consideration. 

6. The outlet water temperature, T^ 2 ' checked to 

see if it exceeds the saturation temperature. 

7. This set of calculations is performed until T ^2 exceeds 
the saturation temperature. This naturally means that 
boiling begins in the last pass of the current total number 
of passes. Initiation of boiling will be treated in the 
boiling section of the model. Therefore, the heating section 
ends with the pass prior to that which initiates boiling. 

The results of the above calculations are: (a) The total 

niimber of passes contained in the heating section, (b) the 
new gas temperature at the inlet to the heating section ' 
and (c) the new water outlet temperature T^ 2 * 

With these new heating section terminal temperatures an 
average outside wall temperature for the section may be 
calculated from 
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wo 



- "fb> 



where total inside/outside areas for the section 

and are the new gas/fluid bulk temperatures for the 

section. The expression for T „ is derived from the follow- 
ing formulation for heat transfer in the heating section. 



Q 





T ^ - T 
gb 



which reduces to 



wo 



= T 



gb 



Ir 



th 






where 



R 



n . h 

t g to 



and 



Ja 



th 



U . . 

oi ti 



The gas-side film temperature is 



T 



gf 



T _ T 
gb wo 
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This gas-side film temperature is now introduced at the 
beginning of the calculations for the heating section, to 
replace the gas bulk temperature. All calculations are 
performed again with 

The gas-side pressure drop in the heating section may 
now be calculated using the previously obtained friction 
factor, f, from the correlation contained in Ref. 7. The 
gas-side pressure drop is 

2f , y - T . 

Ap = max t/r ^ *^w ^ 0 . 14 



A . 
mxn 

= number of tubes per row 
= gas density 

= gas viscosity at wall temperature 
= gas viscosity at bulk temperature. 

D. BOILING SECTION 

The boiling section is solved pass-by-pass since most 
of the section will involve two-phase flow on the cold side 
and the inside heat transfer coefficient will change with 
quality, x, and heat flux, q". The correlation for the 
two-phase region inside heat transfer coefficient selected 



where 



max 



t/r 



w 
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for this model is that recommended by Tong [Ref. 10] for 
both nucleate boiling and forced convection. The equation 
for the two-phase flow heat transfer coefficient was 
given by Schrock and Grossman in Tong [Ref. 10] as 




B[ 



G h 



fg 



+ A( 




n 

) ] . 



The constants are given by Wright [Ref. 10]: B = 6.70x10^, 

-4 

A= 3.5x10 , n = 0.66. The heat transfer coefficient 

assiaming a total liquid flow is 



K 



= 0.023 

d . 

1 



Z r^i ^ 0 . 8 



K, 






and the Martinelli parameter is 



1 



X 



tt 




0.9 




0.5 




0.1 



The limits of the data for this correlation are 



quality; 0.05 - 0.57 

heat flux: 6.0x10^ - 1.45x10^ BTU/hr-ft^. 

In this model, heat flux in the boiling section will be at 

4 2 

levels considerably below 6 x 10 BTU/hr-ft and the full 
range of quality must be considered. Therefore, in order 
to observe the performance of the Schrock and Grossman 
correlation at lower heat flux and quality above 0.57 a 
plot was made (fig. 6) of the ratio h^^^/h^ vs. quality 
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where is the heat transfer coefficient from the Dittus- 
Boelter correlation for water at saturated liquid condi- 
tions. The correlation of Dengler and Addoms as given in 
Tong [Ref. 10] for entirely forced convection vaporization 
is included for comparison. Tong [Ref. 10] indicates that 
low heat flux and high water-vapor mixture velocity favor 
the forced convection mechanism. After a review of the 
information available in Tong, it is not entirely clear 
whether the forced convection mechanism or the nucleate 
boiling mechanism is dominant for the relatively low flow 
velocities and low heat flux of this model. In any case, 
the Schrock and Grossman correlation applied at low heat 
flux represents a type of compromise between the entirely 
forced convection assumptions of Dengler and Addoms and 
the "mixed" assumption of Schrock and Grossman. When the 
Schrock and Grossman correlation is applied under the 
conditions of very low heat flux and quality (x < 0.05), 
as in the first pass of the boiling section of this model, 
it can be seen that the ratio becomes less than 

one^ This is both unrealistic and computationally diffi- 
cult to handle. Therefore, in the boiling section, the 
model assumes that h^p^ = h^ for the first pass or for any 
pass where x < 0.05. This assumption will produce a more 
realistic but still conservative design. 

One other simplification should be explained prior to 
discussion of the actual model calculations for the boiling 
section. Since the boiling section calculations are performed 
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for one complete pass at a time, the average quality, x , in 
a particular pass is used as the "local" quality to compute 
the inside heat transfer coefficient from the Schrock and 
Grossman correlation. 

The first pass in the boiling section, as indicated 
in the calculations for the heating section, will, in 
general, involve both heating water to the saturation tem- 
perature and the initiation of boiling. From the calcula- 
tions for the last pass of the heating section we have the 
following quantities for the first pass of the boiling 
section: water inlet temperature ' average gas 

outlet temperature (T^ ) . The gas side heat transfer 

^out 

coefficient, h^, is obtained, using the average gas tem- 
perature in the boiling section, in the same way as for 
the heating section. 

A rough estimate for the gas temperature into the first 
pass must be obtained in order to begin calculations for 
that pass. This is accomplished by calculating a rough 
overall heat transfer coefficient, for the pass using 

the tube wall resistance, h , and neglecting inside resis- 
tance. A pass effectiveness, is then calculated assuming 

that only boiling takes place in the pass. The expression 
for effectiveness for fluids unmixed between passes for 
C = 0 becomes 



£ 



P 



(1 - exp(-N) 
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where N = NTU- Using e calculated from this expression, 

IT 

the rough gas temperature for the first pass is 






T 



1 

*?out 



Since this pass has been simplified by allowing the inside 

heat transfer coefficient to be computed by the Dittus- 

Boelter correlation this initial guess for is not 

^in 

essential. It would, however, become important if we were 
attempting to "split" the pass with two inside heat transfer 
coefficients, one for heating and one for boiling- Instead, 
the inside heat transfer coefficient, h^, is calculated for 
the entire pass using the Dittus-Boelter correlation with 
water conditions corresponding to the average water condi- 
tions in the heating section of the pass, that is. 



Tf. + T 
^in satf 



■ avg 



With this inside heat transfer coefficient, an expression 
for the overall heat transfer coefficient for the first 
pass may be written. 



U‘ 



oi 






+ 



T: ln(d /d. ) 

ip o 1 

2 u K . L 
w t/p 



A. 



v,l 

op 



The pass effectiveness, e"^, is now calculated in the same 

P 

manner as previously described. With e , the average gas 

IT 
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temperature in is calculated from 



C . 

p min 



- C 

fjn max gput 



in 



C - C 
p min max 



The pass heat transfer may be calculated from 



= m C (T - T ) 

g pg g^_j^ “^out 



The expression for out let enthalpy from the pass is 



= h: 



’out 



in 



_£ 

fh^ 



and outlet quality is calculated from 



h 



out 



^out 



- h 



satf 



fg 



In order to determine the location of the interface between 
heating and boiling the total inside area devoted to heating 
in the first pass is calculated from 



1 , ( !!satf 

hf - hf 
out in 

Subsequent fluid passes in the boiling section are 
calculated until the quality at the outlet of a pass is 
found to be greater than 1.0. When this occurs, the boiling 
section is allowed to end with the last pass for which 
Xout £ A "mixed" pass where both boiling and superheating 
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takes places will be calculated in the superheating section 
of the model. 

After the first pass, boiling section fluid passes are 
calculated in the following manner. 

1. The steam quality at the outlet of the previous pass 

becomes the inlet quality for the current pass, The 

average gas temperature at the inlet to the previous pass 

becomes the outlet average gas temperature, T , for the 

^out 

current pass. 

2. Using the last overall heat transfer coefficient calcu- 
lated for the previous pass, approximate values of 

NTU and e are calculated. 

P 

3. With this approximate the average gas temperature 
into the current pass is calculated from 



T 



in 



T, 




4. Using this gas temperature in, T , the pass heat 

*^in 

transfer, Q , is calculated. This allows the calculation 

of enthalpy out, h^ , and quality out, , as well as 

^out 

the average quality in the current pass, x . 

5. The inside heat transfer coefficient, 
calculated. 



h 



tpf 



Q /A. , 



fg 



tt 



/ 
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where the Martinelli parameter is calculated using 

^tt 

average quality, x , and all other terms are as previously 

3 . 

described. 

6. A new overall heat transfer coefficient is now calcu- 
lated from which a new pass effectiveness can be calculated. 

7. This new pass effectiveness yields a new average gas 
temperature into the pass from which a new pass heat 
transfer may be calculated. 

8. The new pass heat transfer is compared with the pre- 
viously calculated pass heat transfer. If the two heat 
transfer calculations do not agree to within 5% the entire 
set of calculations is repeated. 

9. When the pass heat transfer calculation converges, the 
final quality at the outlet of the pass is calculated. If 
the quality is less than 1.0 another pass is calculated for 
the boiling section. 

10. If the quality is 1.0, the boiling section ends with 
the current pass. If the quality out exceeds 1.0 the boiling 
section ends with the previous pass. 

Calculations for the average gas film temperature and 
gas-side pressure drop are performed in the same manner as 
for the heating section. 

E. SUPERHEATING SECTION 

The first pass of the superheating section will, in 
general, involve both boiling and superheating. The inside 
heat transfer coefficients will be quite different for the 
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boiling and superheating portions for the pass, and, for 
this reason, the pass will be "split" to allow more accurate 
calculation of the heat transfer in these two regions. An 
outline of the calculation scheme is as follows : 

1. Calculate the gas-side heat transfer coefficient for 
the section in the same manner as for the previous two 
sections . 

2. Estimate the average gas temperature into the first 

pass. This is only a rough guess for T , to begin the 

*^in 

calculation. 

3. Calculate the required heat transfer in the boiling 
portion 



Q 



rb 



(h- - h ) 

f f. satv 

in 



where h^ is the enthalpy out of the last pass of the 
^in 

boiling section and is the enthalpy of saturated 

vapor . 

4. Using and the average quality in the boiling portion 

of the first pass, calculate the inside heat transfer 
coefficient, using the Schrock and Grossman correla- 

tion, as in the boiling section of the model. This, together 
with the gas-side heat transfer coefficient and the wall 
resistance, will yield an overall heat transfer coefficient 
for the boiling section, 

5. The calculation of the area required for the boiling 

portion of the pass is performed in an iterative scheme as follows. 
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a. An initial estimate for the required boiling area, 
Al, is made. 

b. The mass flow rate of the gas over the boiling 
portion is then 



• b 

m 

g 



Al m 

A. g 

ip ^ 



c. An effectiveness for the boiling portion may be 
calculated from 



'rb 



C . (T 
mxn '■ g 



- T, 



in 



’sat 



where 



. = m C . 

min g pg 

d. Using the previously calculated overall heat 
transfer coefficient, the estimated boiling area and the 
consequent gas-side heat capacity, NTU is calculated from 



NTU = 



U^. Al 
oi 

=min 



e. A "test” effectiveness, can be found by applying 

this NTU in the relationship for effectiveness 



= 1 - exp (NTU). 



f 

f. The two effectiveness calculations, 
compared . 







are 
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(1) If is greater than then the original 

estimate for the boiling area, Al, was too large for the 
heat transfer required and the specified gas temperature 
into the pass. So, the estimate for the boiling area is 
decreased and a new Al is established. 

(2) If £^ is less than £, , Al is increased. 

b b 

g. A new gas-side heat capacity is calculated from 

ob Al . ^ 

C . = T m C 

min A^p g pg 

h. A new boiling portion effectiveness is calculated 
using the latest boiling area estimate 



^rb 



. (T 
min g 



- T, 



in 



’sat 



i. As before, an NTU is calculated with the new Al. 

This generates a new "test” effectiveness, which is 

compared with £j^. If £j^ and are not equal, a new Al is 
calculated based on the comparison. 

j. These calculations are continued until £, - £^. 

b t 

The current Al is then set equal to A^. 

6. Once the area required for boiling, A^, for a particular 

average gas inlet temperature has been established, the 

following first pass quantities can be calculated; 
sh 

a. A. = area available for superheating 
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b. 



c. 



• b 



m 



g 




d. 



’out 



e. 



T 



sh 

^out 



gas flow rate over the boiling portion 

gas flow rate over the superheating portion 

= gas temperature out of the boiling portion 



= T, 



’out 



g. 

in 



^rb 

C 

g pg 



gas temperature out of the superheating 
portion 



A. T 



T 



sh 

^out 



- 

"^out ^ ^out 

1 



(The weighted average of T 



sh 

^out 



and T 



must equal 



’out 



the T 



previously calculated.) 



’out 



7. Since the heat transfer to the cold side must satisfy 
an energy balance involving the average temperatures of the 
gas, the expected heat transfer in the superheating portion 
may be calculated from the heat release on the gas side 




m 

g 





T 

g 



) 

out 



/ 



and the required heat transfer in the superheating portion 
as 



Q 



sh 
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8. The enthalpy and temperature of the steam out of the 
first pass may be calculated from the assumed heat transfer 



in the superheating portion, 

9. The steam properties for the superheating portion may 

now be found. These properties, along with h^ and wall 

resistance, will yield an overall heat transfer coefficient, 
sh 

for the superheating portion. The Dittus-Boelter 
correlation is applied to calculate the inside heat trans- 
fer coefficient. 

10. The heat capacities on the gas and fluid sides are 

calculated, using the current gas mass flow rate over the 

• sh 

superheating portion, m^ , to calculate the gas heat capacity. 
The heat capacities along with and yield the effec- 

tiveness for the superheating portion of the pass, 

11. A new average gas temperature into the pass may now 
be calculated from the superheating portion effectiveness 
and the superheating portion heat transfer required to 
satisfy the original heat balance from the expression 



= T, 



'sh 



in 



in 



^sh min 



★ 

This temperature, T , is a measure of the temperature 

^in 

which would be required to produce the originally specified 

Q , under the actual conditions of the heat transfer char- 
sh 

acteristics of the superheating portion. 
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12. If the T is the same as the originally specified 
^in 

average gas temperature into the pass, T , then the heat 

*^in 

transfer in the superheating portion, is possible 

under the current distribution of area between boiling and 

superheating and the heat balance is satisfied between steam 

and gas sides for the specified gas temperature in. In 

this event, the calculations for the first pass of the 

superheating section are complete. 

* 



13. If T is less/greater than the current T the 

^in ^in 

average gas temperature into the pass is increased/decreased 

and the entire set of calculations is performed again until 

•k 

T = T . 

g . g . 

^in ^xn 

14. The second and subsequent superheating section passes 
are calculated in much the same way as in the heating section. 
That is, an overall pass effectiveness is calculated for a 
specified number of passes, using the multipass effectiveness 
relationship from Kays and London [Ref. 9] . 

15. The average steam-side heat transfer coefficient, h^, 
is calculated using the Dittus-Boelter correlation with 
steam properties at steam bulk temperature for the remainder 
of the superheating section. The overall heat transfer 
coefficient and the average pass effectiveness, e , are then 

It 

calculated. 

16. The formulation for overall effectiveness with the 
number of passes, n, variable, is 



54 



1 



£ 

oa 





C /C 
min maXxH 



C . /C 



- £ 



min max j n 



C . 
min 

C 

max 



As in the heating section this expression is used in an 
iterative fashion, beginning with n = 1. 

17. Each time e is solved, the gas temperature into the 

oa 

superheating section is calculated from 



T 



g. 

in 



e 

oa 



C . - C T 

min f . max q 

in ^out 

e C - C 
oa min max 



where is the temperature of the steam into the second 

^in 

pass of the superheater and T is the average gas tem- 

^out 

perature out of the second pass of the superheater. This 
gas temperature in is compared with the originally specified 
gas temperature into the heat exchanger, 

18. If the calculated gas temperature in, T , is less 

^in 

than the number of passes is increased by one and the 

calculation is performed again. If T exceeds T ^ , the 

^in ^ 

superheater is allowed to stop with the previous pass. 

Calculations for the average gas film temperature and 
the gas-side pressure drop are performed in the same manner 
as for the heating and boiling sections. 



F. I4ATCHING SUPERHEATER GAS INLET TEMPERATURE 

The model calculations thus far lead to a waste heat 
recovery unit design for which the gas temperature at the 
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inlet is, in general, lower than that specified in the initial 
conditions. This condition occurs because the superheater 
ends with the last complete pass calculated for which the 
entering gas temperature is less than or equal to the 
originally specified WHRU inlet temperature. This outcome 
is, of course, unsatisfactory since the inlet gas temperature 
is an independent quantity which is a function of the horse- 
power setting on the gas turbine. An iterative scheme is 

employed for matching the T calculated in the superheating 

^in 

section with the T^^^ specified in the initial conditions. 

Since the model does not allow the calculation of frac- 
tional passes and since the scheme for the calculation of 
superheater passes will not allow inclusion of a pass which 
would require a gas inlet temperature higher than that 
specified in the initial conditions, the final gas inlet 
temperature calculated for the WHRU will generally be lower 
than that initially specified. In order to match this 
superheater gas inlet temperature with that specified 
initially, the heat exchanger gas outlet temperature is 
allowed to rise. Each time the gas outlet temperature is 
increased, the entire waste heat recovery unit calculation 
is performed again, and the calculated superheater gas 
inlet temperature is compared with the gas turbine exhaust 
temperature. When a match is achieved the design is 
fixed. 
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G. OFF-DESIGN CALCULATIONS 



Once a particular design has been selected the designer 
must test that design at several off-design points. A 
design will define the following WHRU characteristics; 
Dimensions (length, width, height) 

Operating pressure 
Superheater outlet temperature 
Fin-tube configuration. 

The off-design calculation requires that the selected 
WHRU design performance be predicted for a gas inlet tempera- 
ture and flow rate different from that which was used to 
produce the selected design. The design model can be used 
in an iterative fashion to make this calculation. All 
model calculations are made as prescribed in the previous 
sections, with the exception of the pinch point calculations 
which are not performed. 

The off-design procedure requires that the designer 
fix the WHRU frontal dimensions, operating pressure, fin- 
tube configuration and minimum gas outlet temperature. A 
design is then produced for a gas inlet temperature and 
flow rate at an off-design point. This design is checked 
to detemine whether it conforms to the height dimension 
of the WHRU designed at the gas conditions of the off-design 
point. If the height dimension is not matched then the 
number of WHRU passes must be adjusted until a match with 
the design-point WHRU is achieved. The number of WHRU passes 
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may be altered by either adjusting the WHRU gas outlet 
temperature or the superheater outlet steam temperature. 

The effect of either of these adjustments is to change 
the steam flow rate and the heat transfer rate. It would 
be possible to adjust the steam flow rate and allow the 
superheater steam outlet temperature to remain fixed. 

However, for the method of this model, the computation time 
required for this means of control would be large. The 
actual procedure used produces designs relatively quickly 
which closely approximate the results which would be 
obtained by a fine control on steam flow. Each time the 
superheater steam outlet or gas outlet temperature is adjusted, 
a new design is produced. This procedure is continued 
until a design is produced for which the height (number of 
passes) equals that of the design point WHRU. When the WHRU 
height at the off-design gas conditions is equal to that 
of the design point WHRU, the gas inlet temperature is 
matched with the gas inlet temperature specified for the 
off-design point, using the procedure described in Section E. 

H. COGAS SYSTEM OUTPUT MODEL 

The COGAS system output model calculates performance 
parameters for the combined steam and gas turbine system for 
a specified gas turbine input power. The WHRU model pro- 
vides a steam flow rate, pressure, and temperature for a 
specified set of gas turbine exhaust gas conditions which 
correspond to a particular gas turbine horsepower setting. 
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This WHRU output is applied in a simple Rankine cycle to 
calculate the steam-side power output. The gas-side 
pressure drop calculated in the WHRU model is used to 
arrive at a revised gas turbine horsepower. By combining 
these two power outputs, a COGAS system power can be calcu- 
lated. Other performance indicators such as specific fuel 
consumption, thermal efficiency and steam turbine share of 
the load are also calculated in the COGAS system output 
model. 



1 . Rankine Cycle 

The following conditions are assumed for the Rankine 
cycle calculations: 

a. No line losses 

b. Steam turbine efficiency = n = 0.85 

c. Condenser pressure =2.0 psia (4 in Hg) 

d. Feedwater heater pressure = 15 psia 

e. Feedwater temperature = 200 F 

f. Pumping power required is neglected 

g. Fuel LHV = 18400 Btu/lbm 

The following input conditions are received from the 

WHRU model: steam turbine inlet pressure (P^^) , steam turbine 

inlet temperature (T^^) , steam flow rate (m^) , gas turbine 

exhaust pressure drop (AP^.^^) . From P^ and T^^ the turbine 

inlet enthalpy (h^^) and inlet entropy (s^^) can be found. 

The turbine exhaust steam quality, assuming isentropic 

expansion, x- , is calculated from 
^ s 
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fg 
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where 



= entropy of saturated water 
s^g = entropy increment for evaporation. 

The steam turbine exhaust enthalpy assuming isentropic 
expansion, may now be calculated from 



‘2 s 



hf + 



X_ h. 

2s fg 



where 



h^ = enthalpy of saturated water 

h^g, = enthalpy increment for evaporation, 



The isentropic turbine work is 



= h, - h_ , 
ts 1 2s 



and the actual steam turbine work is calculated from 



= n ^ 

t St ts 



The actual turbine exhaust enthalpy and quality are 



*'2 = 4 - ”t 
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and 



X 



2 




In this model it is assumed that heating steam for the 
feed-water heater is provided, via a reducing station, from 
the main steam line. The fraction of the steam mass flow 
rate from the WHRU required for the feedwater heater is 



m = 





where 

h^^ = enthalpy of the feedwater 

= enthalpy of the condensate at feedwater 
heater pressure. 



Therefore, the steam turbine power output is 



P . = (1 - m) m^ W, . 

st ft 

2 . Gas Turbine Performance Calculations 

The original gas turbine horsepower, exhaust gas 
flow rate, and exhaust gas temperature inputs to the WHRU 
model were hand-calculated using figure 4 of Ref. 11 and 
figures 1, 6, and 8 of Ref. 10. The assumptions used to 
calculate the WHRU model inputs were: 
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Ambient temperature = 100 F 

Ambient pressure = 14.696 psia 

Humidity = 0.0 

Fuel LHV = 18400 Btu/lbm 

Inlet loss = 4.0 in. H 2 O 

Exhaust loss = 6.0 in. H 2 O 

In the COGAS system output model, the gas turbine horse- 
power correction factor for the gas-side pressure drop in 
the WHRU (BHP for 6" loss/BHP for WHRU loss) is found by 
solving 



= 1.0125 + 0.002125 AP 

bhp gas 



where = BHP for 6" loss/BHP for WHRU loss. This rela- 

tionship was found by extending the duct loss relationship 
of fig. 7, Ref. 11 linearly from the point of 6 in. H 2 O 
exhaust duct loss already assumed. The final gas turbine 
horsepower is 



gt 






'bhp 



where P^^ is the input gas turbine horsepower. The total 
COGAS system horsepower, is found by adding the steam 

turbine and final gas turbine horsepowers. 

In order to calculate the gas turbine specific fuel 
consumption fig. 4 of Ref. 11 is entered with the final 
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gas turbine horsepower, and the s.f.c. is read from the 

propeller law curve superimposed on that figure. The s.f.c. 

correction factor for exhaust loss in the WHRU, C ^ , is 

sf c 

found in the same manner as C, , from fig. 6 of Ref. 11 

bhp ^ 

where 

= (s.f.c, for total exhaust AP)/(s.f.c. for 6" exhaust AP) 

and the gas turbine specific fuel consumption is 

s.f.c.^, = C - sfc_ . , 

gt sfc 6 in. loss 

The specific fuel consumption for the combined system is 
calculated by first finding the fuel consumption rate of 
the gas turbine 

m- T = s.f.c. ^P,. 
fuel gt gt 

The COGAS system s.f.c. is then calculated from 

_ ^ _ fuel 

'COGAS P. , 

tot 

To allow the direct comparison of the COGAS system with the 

all gas turbine system the specific fuel consumption of the 

★ 

gas turbine at the new, higher COGAS system power, s.f.c.^^, 

is calculated in the same manner as s.f.c. , . 

gt 
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The thermal efficiencies for the gas turbine at 



input power, the COGAS system, and the gas turbine at COGAS 
system power are calculated from the relationship 

2545 

"^th s.f.c. LHV • 



Additionally, the steam turbine share of the load 
is calculated from Pg^/P^ot* 
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III. RESULTS AND CONCLUSIONS 



A. BACKGROUND 

The foregoing model was applied in a computer simulation 
program written in FORTRAN IV for the IBM 360-67 computer. 

A listing of this program appears in Appendix A. Addition- 
ally a set of supporting programs was written for water, 
steam, and air properties. A listing of this set of pro- 
grams appears in Appendix B. 

In order to test the model and develop an understanding 
of the behavior of the WHRU, a set of designs were produced 
for a variety of initial conditions. These designs were 
then reviewed, and an analysis of the model behavior was 
performed. After completion of the model analysis, several 
of the designs were selected for observation at off-design 
conditions. The off-design simulations were produced for 
gas turbine input powers corresponding to ship speeds of 
9, 16, and 20 knots, which represents the entire range of 
gas turbine input powers considered in this thesis. Finally, 
two designs were selected for detailed testing at off-design 
conditions. These latter two designs were tested at gas 
turbine input horsepowers corresponding to speeds of 9 to 
20 knots, at one knot increments. 

B. DESIGN VARIABLES 

The following set of design variables was selected to 
produce the basic set of designs. 
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1. WHRU operating pressure: Three WHRU operating pressures 

of 400, 600, and 800 psia were selected for design production. 

2. WHRU geometric scale: In order to test the model 

response to a range of fin-tube sizes, designs were produced 
for the basic fin-tube configuration described in the geometry 
subsection at 1.0, 0.75, and 0.50 scale. 

3. Frontal dimensions: The approximate WHRU space constraints 

for a DD-963-type engineroom are shown in figures 7, 8 and 

9. Consistent with these constraints, frontal dimensions 
of 12' X 12' and 12' x 15' were selected for the initial 
design set. 

4. Gas turbine input power: The gas turbine input horse- 

power determines the input gas conditions, temperature and 
flow rate, for the WHRU. Since the most advantageous use 
of the COGAS system was assumed to be for "cruise" speed 
conditions, the range of gas turbine input horsepower con- 
sidered corresponds to the gas turbine power required for 
a ship speed range of 9 to 20 knots, with the gas turbine 
operating alone. The three gas turbine input powers selected 
for the production of WHRU designs are shown in Table I. 

Gas Turbine Input Powers 



GTHP 


Approx. Speed (kts) 


Temperature (,F) 


Flow Rate (Ibm/hr) 


16421 


20 


849 


407589 


8526 


16 


742 


328641 


1684 


9 


689 


159731 



Table I 



66 




ETH 9 Rm 

YCNT . 

ry«. C-ANS 



ENGINEROOM 



FIGURE 7 



67 




ENGINEROOM 



FIGURE 8 




SECTION 



ENGINEROOM 



FIGURE 9 



69 



m 



I 



5. WHRU gas outlet temperature: The WHRU gas outlet 

temperature is fixed at 400 F. Reference 3 states that 
the principal cause of "hot" corrosion in the WHRU would 

be due to the formation of H 2 SO^ from an expected concentra- 
tion of SO 3 of up to 50 parts per million. Reference 13 
states that dew point condensation of vapor in the exhaust 
gas becomes a problem for corrosion if the stack temperature 
is depressed below 250 F. Therefore, a minimum average gas 
outlet temperature of 400 F is considered to be reasonably 
safe from the aspect of gas-side corrosion prevention. It 
should be further noted that the outlet gas temperature 
calculated by the program is an average temperature, based 
on the heat release from the gas. Therefore, the actual 
temperature of the gas in the stack will be cooler than 
400 F in some locations due to the non-uniform distribution 
of the gas temperature as it leaves the WHRU. It should 
further be noted that, although this temperature is a fixed 
design variable for design production, it will vary somewhat 
from 400 F in response to the program calculations for 
setting the minimum pinch point and for matching the gas 
inlet temperature with that specified in the initial 
conditions . 

6 . Superheater outlet steam temperature: The target WHRU 

superheater outlet steam temperature is 650 F. This tempera- 
ture will also vary somewhat in response to steam flow rate 
changes brought about by the program calculations to match 
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the gas inlet temperature with that specified in the 
initial conditions. 

7. WHRU water inlet temperature: The WHRU water inlet 

temperature is fixed at 200 F. 

8. The other major COGAS system parameters are fixed as 
follows : 

Condenser pressure: 4 in. Hg 

Steam turbine efficiency: 0.85 

Feedwater heater pressure: 15.0 psia 

Fuel LHV: 18400 Btu/lbm. 

C. THE DESIGN SET 

A set of 54 designs was produced for the combination of 
design variables just discussed. The array of design 
variable combinations is described in Table II. The summary 
computer output pages for the results of each design run 
are presented in Appendix C. 

As discussed in Section A, the purposes for producing 
this design set were to test the model and to promote an 
understanding of the likely behavior of a WHRU for a range 
of conditions. To this end, the following set of output 
variables was considered. 

1. WHRU output 

a. height (number of passes) 

b. gas inlet temperature 

c. gas outlet temperature 
,d. steam flow rate 
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WHRU DESIGN COMBINATIONS AND DESIGN RUN INDEX 
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TABLE II 



e. gas side pressure drop 

f. pinch point temperature difference 

g. heat transfer rate 

2. COGAS system output 

a. gas turbine horsepower (including loss from the 
WHRU gas-side pressure drop) 

b. steam turbine horsepower 

c. COGAS system horsepower 

d. steam turbine share of the load 

e. specific fuel consumption 

(1) s.f.c. of COGAS system 

(2) s.f.c. of gas turbine at COGAS system horsepower 

f. thermal efficiency 

(1) '^t±i COGAS system 

(2) '^th turbine at COGAS system horsepower 

Clearly, several of these variables are related, such as 
s.f.c. and and no additional information is provided by 

studying all variables in a closely related group. Therefore, 
the following subset of output variables was selected for 
analysis . 

1. WHRU height 

2. Steam flow rate 

3. Gas-side pressure drop 

4. Pinch point temperature difference 

5. WHRU heat transfer rate 

6. WHRU gas outlet temperature 
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7. Steam turbine horsepower 

8. COGAS system horsepower 

9. COGAS specific fuel consumption 



The above output variables are presented in tabular form 
using the format of Table II. Table II also serves as an 
index for the set of basic designs. Additionally, the 
significant trends in the variables are presented graph- 
ically. Where the trend of the results was ascertained to 
be the same for the two frontal dimension sets, only the 
values for the 12' x 12' front were presented graphically. 
Where clarity was better served, only one scale was plotted 
as representative of the trend described. 

Effect of Scaling . 

The only significant trends noted with respect to scaling 
were a decrease in WHRU height with a decrease in scale 
and the same trend for gas side pressure drop (see Tables 
III and IV) . The trends are shown in figures 10 and 11. 

The basic fin-tube dimensions of the model were scaled 
by factors of 0.75 and 0.50. For any scale with constant 
frontal dimensions, the total inside and outside areas 
remain the same. That is, for a constant length tube, the 
scaled inside/outside area (A^/A^) is proportional to the 
full scale areas (A^/A^) with the scale factor, s, as the 
constant of proportionality. Hence 
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TABLE III 
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TABLE IV 



HEIGHT VS. SCREE 
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FIGURE 10 
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GRS PRESSURE DROP VS. SCREE 
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FIGURE 11 
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since constant frontal dimensions are to be maintained 
the number of tubes per row is scaled up by the scale 
factor, so 





A. 

1 



and 

op 




S A 

o 
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where is the number of tubes per pass. Therefore, the 

scaled outside and inside areas for a pass for constant 
frontal dimensions are the same as the full scale areas. 
Now, the scaled cross-sectional area for fluid flow on the 
inside will be. 



A 



s 

ff 





s A 



ff * 



So the scaled Reynolds number will be 

s , *■£ ^ =^1 ^ *f <^1 

e s p Ajj p 

which is equal to the Reynolds number for the unsealed 
geometry. It can also be shown that the outside Reynolds 
number remains constant with scaling. Recalling that the 
Dittus-Boelter correlation was used to calculate the inside 
heat transfer coefficient for the single phase regions and 
used with a multiplier for the two phase heat transfer 
coefficient and assuming constant fluid properties, the 
Nusselt number remains constant 
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C . 



N = 0.023 R®*® 

u e r 



So, the scaled inside heat transfer coefficient 



h 



s 

f 



K 

s d . 
1 



C 



C 

s 



increases with a decreasing scale factor. Similarly, for 
the scaled gas-side heat transfer coefficient, 




j 




pl/3 



c 



an increase is seen with a decreasing scale factor. 
Therefore, in addition to a decrease in height attributable 
to a scaling down of the fin-tube dimensions, there is also 
an improvement in heat transfer which results in a further 
reduction in height by elimination of passes. This trend 
is supported by the results . 

Recalling the relationship used to calculate the gas- 
side pressure drop, 

2f , y n 1 H 

AT. - niax t/r , w^O-14 

“g = P 

it can be shown that the scaled gas-side pressure drop is 
inversely proportional to the scale factor. 



AP 



s 

g 



c 

s ' 
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i 



1 

1 



j 



where C is a constant based on the assumption of constant 
gas properties between scales. For the same number of 
passes, then, an increase in gas-side pressure drop could 
be expected. However, since the number of passes is decreased 
with decreasing scale, a net decrease in gas-side pressure 
drop was experienced (Fig. 11) . 

Some differences in pinch point with scale were observed 
in Table V, but no consistent trend could be established. 

It was expected that pinch point AT would remain fairly 
constant with scale for fixed pressure and input gas condi- 
tions. That result was apparent to some extent, particularly 
in the range of high gas temperature and low steam pressure. 
This finding is consistent with the manner in which the pinch 
point AT is established. That is, the program fixes an 
initial pinch point AT only when the rough pinch point AT 
is less than 25 F in the initial calculations. At the lower 
pressures and higher gas temperatures this mechanism would 
not be operative. Also, the final and initial pinch point 
AT calculations are based on the average gas temperatures 
at the inlet to the pass where the AT is calculated. There- 
fore, the actual pinch point may be either lower or higher 
than the one calculated, depending on the actual gas temp- 
erature distribution at the inlet to the pass. The pinch 
point AT will be discussed at greater length in the analy- 
sis of the effects of pressure. 
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PINCH POINT AT [F] 
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TABLE 



1 



Effect of Pressure. 



The most important effect of pressure on the model 
performance arises from the calculations which establish 
the minimum pinch point temperature difference. As related 
in the model description, the minimum pinch point AT is set 
at 25 F. The WHRU outlet gas temperature is adjusted in 
the initial model calculations in order to ensure this 
minimum AT. The initial model calculations establish the 
steam flow rate from a simple energy balance involving the 
gas inlet and outlet temperatures (Tg.^,T^ 2 ) • flow 

rate, the water inlet temperature , the required steam 

outlet temperature , and the steam pressure. Tentative 

interim gas temperatures for the heater inlet, 
boiling section inlet, ^^2' also established in this 

calculation set. In order to establish the minimum pinch 
point AT, the model tests the interim heater gas inlet 
temperature against the criterion 




> T 
- f2 



+ 25 



where T ^2 temperature of saturated water at the 

pressure specified. If T^^ is less than T ^2 ■*" then T^^ 
is set equal to T ^2 25 and a new gas outlet temperature 

is calculated as follows. 



g4 



T T - aT , 
g3 gl 

1 - a 
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where 



a 



^f2 ~ l^fl 
^f4 ” ^fl 



As an aid to understanding the behavior of the model , it 
is useful to predict the WHRU gas inlet temperature at 
which this pinch point calculation becomes operative for 
a particular pressure. This is accomplished by establishing 



where is the temperature of saturated water at the 

pressure specified. The ratio a can also be solved for 
each pressure considered, based on a constant water inlet 
temperature of 200 F and a constant superheater outlet 
temperature of 650 F, 




where h^ = enthalpy of saturated water at the pressure under 
consideration. Finally, assviming a minimum WHRU gas outlet 
temperature of 4-0 F the minimum WHRU gas inlet temperature 
at which the pinch point calculation comes into effect may 
be calculated as follows. 



* 



a "critical" heater gas inlet temperature, » for each 
pressure considered. 
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(1 - a)400 



■gi 



•g3 



a 



Table VI summarizes these calculations for each of the pres- 
sures considered in the design set. 



Minimum Gas Inlet Temperature for Pinch Point Calculation 



Pressure 

(psia) 




* 

a T 1 (F) 





400 


470 


.219 


720 


600 


497 


.262 


770 


800 


543 


.299 


878 



Table VI 



Now the pressure at which the pinch point calculation becomes 
operative for each gas turbine input power considered is 
predicted. The pressure for which the pinch point calcula- 
tion begins to adjust gas outlet temperature is listed below 
for each gas turbine input power considered, in Table Vll. 

Pressure At Which Pinch Point Calculation Takes Effect 



* 



Gas Turbine 
Input Power 


V 


Pressure 

(psia) 




16421 


849 


800 


878 


8526 


742 


600 


770 


1684 


689 


400 


720 



Table VII 
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It should be recalled that the outlet gas temperature is 
also adjusted by the iterative technique employed to match 
the model-calculated gas inlet temperature with the gas 
inlet temperature of the initial conditions, once the 
design is established. This is a relatively minor adjust- 
ment, however, and is never greater than 10 F. These adjust- 
ments to WHRU gas outlet temperature are apparent in the 
values of Table VIII and are displayed graphically for the 
12' X 12' frontal dimentions and 0.75 scale in Fig. 12. 

The final calculated pinch point for each of the design 
sets reflected the influence of the initial pinch point 
calculation. These results were shown in Table V and Figure 
13. It is re-emphasized here that the calculated pinch 
point AT is useful only as an indicator of what the actual 
pinch point might be for a particular heat exchanger design. 
This is true because the model neglects the possibility of 
a gas temperature distribution across the inlet of a V7HRU 
pass and, instead, uses the average gas temperature for 
calculations . 

There are two other trends in the design results which 
emanated directly from the control exercised by the model 
over the pinch point AT and the resulting gas-fluid tempera- 
ture distributions in the designs produced. These trends 
were apparent when the effects of pressure on WHRU height 
and WHRU heat transfer/steam flow rate were examined. 
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WHRU GAS OUTLET TEMPERATURE [F] 
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TABLE VIII 
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GnS OUTLET TEMP. VS. PRESSURE 
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FIGURE 12 



PINCH POINT VS. PRESSURE 
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FIGURE 13 



The effects of pressure on WHRU heat transfer and steam 
flow rate are shown in Figures 14 and 15 respectively. 

Tables IX and X list the actual values obtained for WHRU 
heat transfer and steam flow rate for the designs produced. 
WHRU heat transfer and steam flow rate directly reflected 
the adjustment of the WHRU gas outlet temperature to main- 
tain the minimxim pinch point AT. That is, as the WHRU 
outlet gas temperature was raised, the heat transfer rate 
and steam flow decreased. The small increase shown in 
Figures 14 and 15 in heat transfer and steam flow rate 
for 16421 gas turbine input power at the 0.50 scale reflects 



only a minor adjustment difference in T for matching 

'^out 

T between 400 psia and 600 psia. 

^in 

The relationship of WHRU height to WHRU pressure is 
slightly more complex. Figures 16, 17, and 18 demonstrate 
the approximate gas-fluid temperature distribution along 
the total length of the WHRU for the 0.75 scale at the 
three pressures and input horsepowers considered. The 
general trend observed in these diagrams of decreasing gas- 
fluid temperature difference with increasing pressure is 
representative of all the designs produced. This decrease 
in gas-fluid temperature difference, taken by itself, would 
produce larger heat exchangers with increasing pressure. 
This trend may be observed for the 12 ’ x 12 ' front at the 
high gas turbine input horsepower in Figure 19. At the 
medium and low gas turbine horsepowers, however, the WHRU 
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FIGURE 14 
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FIGURE 15 



HEAT TRANSFER RATE [Btu/hr] 





o 


(j\ 


in 


00 




O 




n* 


VO 


in 




in 


<y\ 




o 


o 


00 




00 




H 




• 


o^ 


VO 


a\ 


VO 


in 


a\ 


CT^ 


H 


cn 




o 


1— ( 


H 


00 


CN 


<Ti 


o 




VO 


VO 




II 


CN 




H 


00 


o 


VO 


(Ti 


o\ 


CN 




<D 


CN 




in 




o 


H 


VO 


cn 


H 


- 


rH 


in 


in 






in 


cn 


o 


a\ 


00 


ID 


ClJ 








CN 


CN 


CN 


H 






tH 


o 






















cn 




















X 






















_ 


m 




















OJ 


(>« 


o 


r- 




cn 


(j\ 


G\ 


H 


H 








rH 


cn 


CN 


cn 






in 


VO 


cn 




o 


00 


in 


a\ 


H 


in 


H 


CN 


<y\ 


CN 




II 


rH 


o 




a\ 


00 


in 


cn 


cn 




• • 


<D 


H 


00 


o 


m 


CTv 


cn 


cn 


CN 




4-> 


fH 


CN 


VO 


cn 


o^ 


H 


VO 


H 


in 


H 


c 


ClJ 


in 


in 




r- 


in 


CN 


H 


as 


00 


o 


o 








CN 


CN 


CN 


H 








CO 




















Cl, 
























o 


o 


(Ti 




(J\ 


o 


m 




CN 


H 






cn 


<y\ 


cy\ 


00 


c\ 


00 




VO 








00 


a\ 


H 


H 


H 


VO 


CN 


o 


cn 




(1 


00 


H 


VO 


in 










o 






00 


CN 


cn 


00 


CN 


o 


cn 


VO 


H 




'«V 


in 


CN 


00 


00 


in 


00 


H 




H 




rrt 


m 


in • 


cn 


r' 


m 


CN 


H 


as 


00 




tu 

o 








CN 


CN 


CN 


H 








CO 






















c 

1 o 


CN 


m 


CN 


cn 


H 


H 




cn 


O 




U J 


in 


a\ 


VO 


r- 


O 


VO 


VO 




O 






VO 


VO 




G\ 


• <J\ 


H* 


CN 


in 


CN 




CJ 


cn 


00 


cn 


VO 


cn 


CN 


in 


o 






II 




H 


VO 






VO 


CN 


VO 


cn 




CD 


r' 


H 




00 


H 




H 


H 


H 




1— 1 




VO 






in 


CN 


H 


as 


00 


CM 


ClJ 

(j 








CN 


CN 


CN 


H 








CO, 




















X 
























in 




















- 


c^ 


CN 


VO 


m 


VO 


<y\ 


H 


VO 


00 


VO 


CM 


• 


a\ 


cn 


a\ 


H 


o 


CN 




cn 






o 


CN 


VO 


VO 


in 










H 




11 


o 




o 


in 


o 


a\ 


r' 


o 


VO 




0) 


cn 


00 






00 


VO 


H 




CN 


• • 


t — i 


H 


a\ 


in 




CN 


in 


O 


CN 


CN 


-P 


cd 


in 




cn 


r' 


in 


CN 


H 


as 


00 


c 


o 








CN 


CN 


CN 


H 






o 


CO 




















Pi 






















Cjl< 
























o 


00 


cn 


VO 


CN 


<J\ 


H 


CN 




H 




• 


cn 




o 


cn 


O 


O 


in 


in 








H 


cn 


rr 




H 


<J\ 


cn 


as 


CN 




II 


cn 


a\ 


cn 


cn 


VO 


in 


cn 




O 




(D 


00 




CN 


m 


00 


H 




CN 






rH 


CN 


o 


CN 






00 


o 


in 


H 




cd 


in 


in 




r' 


in 


CN 


H 


cn 


00 




o 








CN 


CN 


CN 


H 








CO 






















(D 






















P 






















D 


o 


O 


o 


O 


O 


O 


O 


O 


O 




CO 


o 


o 


o 


o 


O 


o 


o 


o 


O 




cn 




CO 


CO 


N" 


CO 


CO 




(O 


00 




CD 






















P 






















PU 


























H 






CO 






N- 






PU 




CM 






CM 






CO 






pp 










in 






CO 










CO 






CO 






H 






o 




H 

















93 



TABLE IX 




4 



STEAM FLOW RATE [Ibm/hr] 
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El CURE 16 



TEMPERATURE DISTRIBUTION DIAGRAMS 
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FIGURE 17 
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FIGURE 18 



HEIGHT VS. PRESSURE 
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FIGURE 



height remained fairly constant or actually decreased with 
increasing pressure. This trend is explained by the 
decreasing heat transfer rate with increasing pressure 
attributable to the pinch point calculation adjustment to 
WHRU gas outlet temperature shown in Figure 12 . 

The design set results for steam turbine horsepower 
output are given in Table XI. These results are presented 
graphically for the 12' x 12 ' front in Figure 20. Two 
separate trends were identified. The steam turbine horse- 
power demonstrated an overall increase with increasing 
pressure for the high gas turbine input horsepower. Although 
there is a small net decrease in WHRU heat transfer rate 
with increasing pressure at the high gas turbine input 
power, the steam is rejected to the condenser at a lower 
quality for the higher pressures. Therefore, there is less 
heat rejection in the condenser for higher pressures. This 
better performance of the high pressure systems in the 
Rankine cycle more than compensates for the small decrease 
in WHRU output at these pressures. For the medium and low 
gas turbine input powers the trend of decreasing WHRU output 
with increasing pressure (Fig. 20) was dominant in the 
Rankine cycle also, producing systems of decreasing steam 
turbine horsepower with increasing pressure . 

The design set results for the COGAS system combined 
horsepower are given in Table XII. The results are displayed 
graphically in Figure 21 for the 12' x 12' front. It is 
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TABLE XI 



HORSEPOWER VS. PRESSURE 
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FIGURE 



COGAS SYSTEM HORSEPOWER 
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TABLE XII 



COGnS HORSEPOWER VS. PRESSURE 
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FIGURE 21 



apparent that the system horsepower output trend followed 
that for steam turbine horsepower output. This was the 
expected result. This trend is modified only slightly by 
additional gas turbine horsepower losses due to a changing 
gas-side pressure drop across the WHRU. 

Finally, the system specific fuel consumption followed 
closely the results for steam turbine horsepower output. 

That is, slightly better systems were observed at lower 
pressures for the medium and low gas turbine input powers 
and a slightly better system was produced at higher pressures 
for the high gas turbine input horsepower. These trends 
are shown in Figure 22 and Table XIII. 

Effect of Gas Turbine Input Power . 

An increase in WHRU height with increasing gas turbine 
input power was noted in the results for the design set 
produced (see Fig. 23 and Table III) . This increase in 
WHRU height generally followed the increase in WHRU heat 
transfer rate with increasing gas turbine input power. Fig. 
24. This trend is modified only by a slightly increasing 
outside heat transfer coefficient with increasing gas turbine 
input power (see calculations for gas-side heat transfer 
coefficient in the model description) . 

As expected, the COGAS system specific fuel consumption 
improved (Fig. 25) with increasing gas turbine input power. 
This improvement represents the trend of improvement in 
the gas turbine itself with increasing power, since the 
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COGRS S. F. C. VS. PRESSURE 
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PRESSURE (PSIR) 
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TABLE XIII 



HEIGHT VS. GRS TURBINE INPUT 
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GRS TURBINE INPUT PONER 



HERT TRRNSFER RRTE VS. POWER 
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GRS TURBINE HORSEPONER 



COGRS S. F. C. VS. POWER 
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GRS TURBINE HORSEPONER 



contribution of the steam system decreases with increasing 
power . 

Effect of Frontal Dimensions . 

The only significant trends for changing frontal dimen- 
sions were in WHRU height and volume (see Tables III and 
XIV) . The expected result of decreased height with increased 
frontal area was observed. The gas-side pressure drop 
follows this decrease in height. The WHRU volume showed a 
significant increase for increased frontal area for all 
designs produced. This result is consistent with the model 
formulation. That is, as the frontal area is increased 
and the tube length is held constant the area for fluid flow 
is increased with a consequent decrease in Reynolds number 
and inside heat transfer coefficient. The same result is 
seen in the outside heat transfer coefficient. As the 
frontal area is increased, the minimum flow area for the 
gas is increased, causing a drop in the gas side Reynolds 
number and heat transfer coefficient. 

D. OFF-DESIGN RESULTS 

The off-design results were produced using the same 
computer simulation program used to produce the design set. 
The method used for off-design runs was described in Section 
F of the model description. As the computer program is 
presently designed, the designer must manually adjust the 
inputs to produce an off-design run. The procedure is 
reviewed briefly below. 
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WHRU VOLUME [ft 
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TABLE XIV 



First, the designer selects a WHRU design to be inves- 
tigated at off-design conditions. Once this design has 
been selected, the pressure, frontal dimensions, height and 
scale of the WHRU are fixed. The designer then enters the 
gas flow rate and gas inlet temperature corresponding to 
the off-design point to be investigated. A design is pro- 
duced with these gas conditions and the physical charac- 
teristics of the WHRU design under consideration. The 
designer then checks the number of passes (height) for the 
resulting design. If the number of passes does not match 
that of the selected design, the WHRU gas outlet temperature 
and/or superheater steam outlet temperature are adjusted 
to increase or decrease the number of passes until a match 
is achieved. Once the physical characteristics of the off- 
design point WHRU match those of the design-point WHRU, the 
gas inlet temperature is matched with the initial conditions, 
and the performance of the selected WHRU design is established 
for the off-design point under consideration. 

In order to test the feasibility of the procedure described 
above and to demonstrate the off-design performance of a 
wide variety of WHRU designs, ten WHRU designs were selected 
for off-design runs. These designs are listed in Table XVI. 
None of the full-scale designs was considered because the 
heights of the designs produced at that scale were, for the 
most part, either too large or borderline according to the 
space limitation described in Figure 9. 
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The results of these off-design runs are given in 
Table XVII. The summary output pages for each of these 
off-design runs along with the design-point runs are pro- 
vided in Appendix D. In terms of specific fuel consumption, 
the performance difference between the design sets con- 
sidered was not large. The performance range for fuel 
consumption is even less significant. Table XV demonstrates 
the performance range for the designs considered. 

WHRU Off-Design Performance Ranges 



Power 


s.f .c 




fuel consump- 
tion (Ibm/hr) 


approx, fuel consump- 
tion difference (gal/hr) 


high 


.329 - 


.340 


6992.8 - 


- 7019.3 


3.9 


medium 


.391 - 


.403 


4429.2 - 


• 4439.7 


1.5 


low 


.615 - 


. 658 


1764.6 - 


- 1767.5 


0.4 



Table XV 

Since the fuel consumption comparison showed very little 
difference between the ten designs considered and since all 
designs considered satisfied the dimensional constraints 
of figures 7, 8, and 9, only two designs were selected for 
additional off-design runs. 

For these additional off-design runs, a pressure of 
600 psia was selected. This pressure did not produce the 
smallest designs but it is a steam pressure with which the 
Navy has some experience in terms of materials requirements, 
maintenance, and operation. A scale of 0.75 was selected 
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WHRU DESIGNS CONSIDERED 



Design 

Run # Pressure Scale 



5 


800 


.75 


14 


600 


.75 


16 


600 


.50 


22 


400 


.75 


24 


400 


.75 


35 


800 


.50 


40 


600 


.75 


41 


600 


.75 


53 


400 


.50 



Design 



Front 


Power 


Height 


12x12 


med 


2.9 


12x12 


med 


3.2 


12x12 


high 


2.0 


12x12 


high 


2.9 


12x12 


low 


2.4 


12x15 


med 


1.5 


12x15 


high 


3.4 


12x15 


med 


2.9 


12x15 


med 


1.5 



TABLE XVI 
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WHRU OFF-DESIGN RESULTS 
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TABLE XVII 



which yields a tube outside diameter of 1.5 inches and a 
fin spacing of about 8 fins/inch. The rationale for this 
selection was that a scale of 0.50 with 1.0 inch OD tubes 
and a fin spacing of 12 fins/inch would probably be too 
susceptible to both inside and outside fouling. In order 
to achieve the closest possible comparison, the design-point 
for both designs selected was for the gas conditions corres- 
ponding to the medium (8526 BHP) gas turbine input horse- 
power. Both frontal areas were used for these detailed 
off-design runs. Table XVIII summarizes the characteristics 
of the two designs selected. 

Final Design Characteristics 
design GT input 



run # 


design point 


scale 


front 


pressure 


height 


14 


8526 


0.75 


12 X 12 


600 


3.2 


41 


8526 


0.75 


12 X 15 


600 


2.9 



Table XVIII 

Table XIX summarizes the performance results for the 
COGAS system with the 12 ' x 12 ' front WHRU and Table XX 
gives the results for the system with the 12 ' x 15 ' front 
WHRU. In both cases, the performance characteristics of 
the gas turbine alone, at the COGAS system horsepower, are 
provided for comparison. The summary output pages for each 
of these runs are provided in Appendix E. 
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PERFORMANCE OF 12'xl2' FRONT WHRU 
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TABLE XIX 



PERFORMANCE OF 12'xl5* FRONT WHRU 
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TABLE XX 



The specific fuel consumption curves vs. brake horse- 
power for the 12' x 12' front and the 12' x 15' front are 
plotted in figures 26 and 27 respectively. The difference 
in s.f.c. performance between the two WHRU designs con- 
sidered was not significant. Geometric considerations were 
therefore used in selecting the better of these two designs. 
A geometric comparison of the two WHRU designs is given in 
Table XXI. 

Geometric Comparison of the Final Two Designs 



Front 


Frontal 

Area 


Height 


Vo lume 


Total Out- 
side Area 


12' xl5' 


181.9 


2. 9 


528.96 


49586.4 


12' xl2' 


144.8 


3.2 


464.64 


42776.5 



Table XXI 



It was concluded that, since both designs met the height 
constraint of figure 9, the 12' x 12' front was the 
better choice because of the significantly lower volume and 
total outside area requirements. 

After the selection of the 12' x 12' front WHRU as the 
final design, estimates were made of the possible COGAS fuel 
savings for a DD-963-type destroyer over 1000 hours of 
operation. The NAVSEC Standard Destroyer Profile (Fig. 28) 
was used to determine the number of hours a destroyer-type 
ship spends at each speed during 1000 hours operation. 
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COGHS S.F.C. VS. BHP (12 'XI 2' FRONT) 
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FIGURE 26 



COGRS S.F.C. VS. BMP (12^X15^ FRONT) 
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FIGURE 27 
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SPEED (KNOTS) 



Figure 29 was used to obtain the shaft horsepower require- 
ments for each speed. Finally, figure 30 was used to obtain 
the s.f.c. of the gas turbine operating alone. 

Two estimates of fuel savings were made. The first 
estimate assumed that the COGAS system would be operated for 
speeds of 5 to 20 knots and would be secured outside that 
range. The second estimate employed the COGAS system for 
speeds of 5 to 23 knots. For the latter operating range, 
the criterion for the high COGAS speed was an imposed maxi- 
mum gas turbine horsepower input of 20000 BHP. Both fuel 
savings estimates are based on the following assumptions 
and simplifications: 

1. The COGAS mode of operation implies the use of one 
engine (COGAS mode) on one shaft with the other shaft 
dragging. 

2. The gas turbine is assumed to operate at idle (.1000 BHP) 
at speeds below 8 knots. 

3. Maneuvering combinations, when two main engines would 
normally be on the line, are not considered. 

4. In the speed range of 21 to 27 knots, in the pure gas 
turbine mode, one engine per shaft is used, with two shafts 
on the line . 

5. From 28 to 32 knots four gas turbines are on the line. 
Table XXII provides the estimates of s.f.c. and fuel 

consumption for COGAS between 5 and 20 knots and for pure 
gas turbine for all speeds. The total estimated fuel savings 
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SHIP SPEED (KTS) 



GflS TURBINE 
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FIGURE 30 
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Hours 


BHP 


COGAS 

s . f . c . 


GT 

s . f . c . 


COGAS 

CONS 


GT 

CONS 


10 










17510 


12.1 










21187 


12.1 










21187 


12.1 


GT 


est. avg. 


est. avg. 


21187 




IDLE 


.800 


1.751 






12.1 


at 






9680 


21187 




1000 










12.1 








9680 


21187 


12.1 








9680 


21137 


12.1 








9680 


21187 


12.1 


1684 


.750 


1.100 


15282 


22414 


12.1 


1895 


.705 


1.070 


16165 


24534 


63 


2105 


.680 


1.000 


90178 


132615 
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3158 
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130760 


173500 


70 


4316 


.560 


. 700 


169187 


211484 


71 


5579 


.512 
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202808 


241626 


68 


7368 


.460 


.550 


230471 


275563 


61 
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232497 
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66 


10421 


.400 


. 489 


275114 


336327 


72 


12316 


.380 
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399372 


481435 


89 


16526 


. 350 


. 433 


514785 


636862 



TABLE XXII 
Page 1 
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22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 



COGAS 

CONS 



GT 

CONS 



COGAS GT 

Hours BHP s.f.c. s.f.c. 



75 


17368 


.519 


676049 


43 


20000 


.495 


425700 


15 


22105 


.482 


159819 


13 


25263 


.463 


152058 


7 


28684 


.448 


89953 


5 


32631 


.435 


70972 


4 


37895 


.418 


63360 


3 


44210 


.482 


62928 


3 


51580 


.460 


71180 


3 


58947 


.444 


78517 


3 


69474 


.427 


88996 


3 
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.415 
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4774897 5429976 



12.1% savings 
over 1000 hrs operation 
655088 Ibm savings 
or 96336 gal. dist. fuel 

TABLE XXII 
Page 2 
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over 1000 hours of operation was 12% which is roughly 
equivalent to 96000 gallons of distillate fuel. 

The estimates of s.f.c. and fuel consumption for pure 
gas turbine and for COGAS for the 5 to 23 knot speed range 
are given in Table XXIII. The total estimated 1000-hour 
fuel savings for this COGAS operating range was 19% which 
translates to approximately 152000 gallons of distillate 
fuel. 

E. CONCLUSIONS 

The waste heat recovery unit model and the accompanying 
simple COGAS system output model provide the basic frame- 
work for additional studies of the application of the COGAS 
system to U. S. Navy gas turbine powered ships. The model, 
as presently formulated, provides a reasonable estimate of 
WHRU size required for the fin-tube configuration considered. 

Estimates of the COGAS system performance in the 5 to 
20 knot range and the 5 to 23 knot range are encouraging. 

Even considering the assumptions involved in making the 
fuel savings estimates, it is likely that a COGAS version 
of a DD-963-type ship would consume on the order of 300000 
to 450000 gallons less fuel in one year of operation than 
the current DD-963 class ships. 

The most probable COGAS system would be one which would 
be designed to operate in non-maneuvering situations for 
either speed range considered above. That is, a system which 
would be employed for maximum fuel savings in the cruise mode. 
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19.0% savings 

fuel saved in 1000 hrs operation 
1033115 Ibm 

or 151928 ga. dist. fuel 
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During maneuvering and high speed operations the gas 
turbines would be operated alone with the waste heat 
recovery unit dry. Because of the relatively low gas 
turbine exhaust temperatures involved, it is not expected 
that damage would occur in the WHRU when operated dry. 
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IV. RECOMMENDATIONS FOR FURTHER RESEARCH 



A. WHRU DESIGN MODEL IMPROVEMENT AND EXPANSION 

The type of control imposed on the WHRU design model of 
this thesis was essentially the adjustment of the steam flow 
rate to maintain a superheater steam outlet temperature of 
650 F at the pressure specified. This steam flow rate 
control was further extended to accommodate the maintenance 
of a minimum pinch point AT of 25 F. Further investigation 
of WHRU design should include the effects of allowing the 
superheater outlet steam temperature and/or pressure to 
"float" with changing gas turbine exhaust conditions. The 
combination of controls on the superheater outlet pressure, 
temperature, and flow rate could lead to improved WHRU output 
while still maintaining an acceptable minimum pinch point 
temperature difference. 

In order to increase the precision of. the model and lead 
to increased knowledge of the gas temperature distribution 
in the WHRU, the present model could be modified for a "fine 
mesh" approach. That is, instead of solving one pass at a 
time, each pass could be divided into a number of segments. 
The solution could then proceed using either the finite 
difference or finite element method. 

Other possible model improvements are as follows: 

1. Include the calculation of fluid/steam-side pressure 
drop. 
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2. Additional fin-tube configurations should be considered 
for possible improvement of the heat transfer characteristics 
and the weight and space requirements of the WHRU. 

3. Consideration of inside and outside heat transfer sur- 
face fouling should be given when investigating possible 
fin-tube configurations. 

4. Working fluids other than water should be investigated 
for possible enhanced thermodynamic characteristics at a 
minimum cost in system maintainability. 

B. COGAS MODEL IMPROVEMENTS 

After the waste heat recovery unit, the steam turbine 
and condenser should have the highest priority for modeling. 
Both of these units were modeled as "black boxes" in this 
model. The modeling or performance mapping of "state of the 
art" steam turbines for application in the COGAS model would 
allow the designer to further refine the estimate of system 
performance for various WHRU outputs of steam temperature, 
pressure, and flow rate. The inclusion of a steam turbine 
model would also allow the designer to make size and weight 
estimates for that component. The same benefits would be 
derived from the inclusion of a condenser model in the COGAS 
system. As the condenser pressure is reduced, the system 
performance should improve. The improvement must, however, 
be weighed against the increased condenser size necessary 
at the reduced pressure for constant steam and cooling water 
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conditions. A condenser model would allow the designer 
to make these comparisons. 

In the model of this thesis, pumping power for the 
condensate and feedwater was assumed to be negligible. In 
the case of the condensate pump, this is probably not a bad 
assumption since it is likely that electric pumps could be 
used, and the cost to the ship's electrical output would 
not be large. The feedwater pump, however, would probably 
be steam turbine driven at some cost to the COGAS system 
output. The required feed pump pumping power for various 
loads should be mapped and included in the model. 

Finally, the likely engineroom machinery layout for the 
COGAS system should be described in enough detail that the 
lengths of piping runs, number of turns, and location of 
valves could be estimated. The inclusion of this feature 
in the model would allow the designer to address pressure 
drop and heat losses between components. This feature would 
also contribute to the COGAS system weight and space predic- 
tion. 

C. SYSTEM OPTIMIZATION 

Since the COGAS system is operated with a relatively 
small heat source, optimization should be attempted in the 
design of the system. The technique of non-linear programming 
could be applied to a COGAS system optimization where an 
important system output parameter, say specific fuel con- 
sumption, is expressed as a function of the system design 
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variables such as component weight, component dimensions, 
and gas turbine backpressure. This function is called the 
objective function, and it is maximized or minimized subject 
to constraints which are also expressed as functions of the 
design variables. These constraint functions may be either 
linear or non-linear and express some parameter limit not 
to be exceeded, such as WHRU or condenser total volume. 

For a system of the complexity of the COGAS system a 
"local" optimization approach would probably have to be 
adopted. That is, each of the major sub-models (WHRU, 
steam turbine, and condenser) could be optimized using local 
constraints. Then the entire system could be optimized using 
linking variables, such as enthalpy, flow rates, and tempera- 
tures between components. The process would be repeated 
as many times as necessary to achieve the final system 
optimization. One available computer-based system using 
the technique of non-linear programming is the COPES/CONMIN 
program [Ref. 14 J . 
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APPENDIX B 



This appendix contains the program listings for the 
subroutines called by the WHRU design program and the COGAS 
model to obtain steam, water, and air properties. Refer- 
ences 15 and 16 provide relationships for use in calcula- 
ting the thermodynamic properties of steam and water. The 
relationships of Ref. 15 for the thermodynamic properties 
of steam were used in Subroutine SS of this appendix. 

The remainder of the subroutines for steam, water and 
air properties were obtained by fitting polynomials to the 
values for those properties found in References 17 and 18. 

In all cases, the polynomial coefficients were obtained 
using a "canned" program for least squares regression pro- 
vided with the Hewlett-Packard System 9845 Computer installed 
in the Naval Postgraduate School Mechanical Engineering 
Department. 
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brake EOHSEPnwERj E526.0, APPROXIMATE C URRE SPONOI ^G SHIP SPEEDS 16.0 KTS 
EXhAlSI CAS TEMPERATURE: 742.0 E 

EXhACSr CAS iLOk RATE: 32B641.0 LdP/HR I 91.3 LEH/SECl 
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GAS URBINe 

RR4Ke HORSEPOWERS U421.0, APPROXIMATE CUKRESPONOUG SHIP SPEEDS 20.0 KTS 
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EXHAIST GAS FLOW RATES 407589*0 L8P/HP U13*2 lEP/SECI 
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